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NOTATION 


Small  cross-sectional  area  of  a  water  column 
Shock  wave  velocity  (•  5000  fps  for  sea  water) 

Base  of  natural  logarithms  (»  2.71828) 

Deceleration  due  to  gravity  («  32.2  ft/sec2) 

N'ondimensional  deceleration  due  to  gravity 

Depth  below  surface 

Depth  below  surface  ( h  -  h  ‘  ^  0) 

Exponential  depth  decay  constant 

Pressure  of  incident  shock  wave  at  a  given  time  after  arrival 
Pressure  in  the  slab  just  above  the  cavitated  zone 
Air  pressure  («.  14. 7  psi) 

Peak  pressure  of  shock  wave 

Time  after  arrival  of  the  incident  shock  wave 

Time  at  which  the  surface  slab  has  grown  to  a  given  thickness 

Time  after  kickoff  (/  ^  t  *  ^  0) 

Closure  time 
Time  at  which  y  =  « 

N'ondimensional  time  at  which  the  surface  slab  has  grown  to  a  given 
thickness 

Nondimensional  closure  time 

Kickoff  velocity  of  the  water  at  a  given  depth 

Kickoff  velocity  of  the  water  surface 

Velocity  of  the  water  surface  at  a  given  time 

Velocity  of  the  water  surface  at  a  given  time  without  the  effect  of 
gravity 

Nondimensional  velocity  of  the  water  surface  at  a  given  time 

Nondimensional  velocity  of  the  water  surface  at  a  given  time  without 
the  effect  of  gravity 

Thickness  of  the  surface  slab  at  a  given  time 
Depth  below  surface  (y  -  y '  -  0) 

Depth  of  water  particle  which  falls  back  to  its  initial  position  before 
joining  the  surface  slab 
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y e  The  closure  depth 

y*  Nondimensional  thickness  of  the  surface  slab  at  a  given  time 

y *  Nond  mensional  closure  depth 

y  j  Depth  of  upper  bound  of  the  below  surface  slab 

y2  Depth  of  lower  bound  of  the  below  surface  slab 

2  Depth  below  surface 

fi  Angle  which  a  plane  wave  makes  with  the  surface 

S  Surface  displacement 

Sg  Surface  displacement  w  ithout  the  effect  of  gravity 

5*  Nondimensional  surface  displacement 

Nondimensional  surface  displacement  without  the  effect  of  gravity 
p  Density  (  *2  slugs 'ft3  for  sen  water) 

r  Surface  cutoff  time  at  a  given  depth 

0  Exponential  time  delay  constant 
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ABSTRACT 


A  theory  of  bulk  cavitation  is  presented  in  which  an  equation  governing 
the  motion  of  the  water  surface  is  derived.  This  equation  is  found  to  have  a 
simple  exact  solution  with  which  calculations  can  be  performed  without  a 
computer.  This  solution  not  only  yields  the  results  of  Walker  and  Gordon  for 
the  closure  depth  and  time,  but  also  predicts  the  complete  surface  velocity 
history  rather  than  just  a  straight  line  model  of  the  surface  velocity  history. 
Simple  exact  equations  for  the  surface  slab  thickness  history  and  the  sur- 
face  displacement  history  are  also  derived.  The  theory  is  found  to  be 
physically  plausible  and  is  in  reasonable  agreement  with  the  experimental 
data  of  Walker  and  Gordon. 

ADMINISTRATIVE  INFORMATION 

The  work  of  this  report  was  funded  under  Defense  Atomic  Support  Agency  Subtask 
NB004(3089). 


INTRODUCTION 

When  a  shock  wave  caused  by  an  underwater  explosion  reaches  the  surface,  it  is  re¬ 
flected  as  a  negative  pressure  wave  which  is  superposed  on  the  incident  positive  pressure 
wave.  At  a  certain  depth  below  the  surface,  the  total  pressure  may  become  negative.  At 
about  this  depth,  the  water  cavitates  since  under  the  usual  test  conditions,  water  seems  un¬ 
able  to  withstand  substantial  negative  pressures,1  The  cavitation  may  persist  down  to  sig¬ 
nificant  depths  and  cause  what  is  known  as  "bulk  cavitation."  A  knowledge  of  this  phe¬ 
nomenon  is  necessary  to  determine  the  motion  of  the  water  and  thus  the  effect  of  underwater 
explosions  on  surface  ships  and  submarines. 

This  report  presents  an  incompressible  theory  of  bulk  cavitation  in  which  an  equation 
governing  the  motion  of  the  water  surface  is  derived.  This  equation  is  found  to  have  a  simple 
exact  solution  with  which  calculations  can  be  performed  without  a  computer.  This  solution 
not  only  yields  the  results  of  Walker  and  Gordon1  for  closure  depth  and  time,  but  also  pre¬ 
dicts  the  complete  surface  velocity  history  rather  than  just  a  straight  line  model  of  the  sur¬ 
face  velocity  history.  Other  theoretical  details  have  been  presented  in  reports  by  Cushing 
et  al.2*  3  A  more  involved  theory  which  includes  the  compressibility  of  the  water  is  being 
studied  and  will  be  the  subject  of  later  work.  |t  is  expected  that  the  results  of  the  com¬ 
pressible  theory  can  also  be  applied  without  the  aid  of  a  computer. 
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Reference*  are  lifted  on  page  11. 
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THEORY  OF  BULK  CAVITATION 

After  a  shock  wave  reaches  the  surface,  ii  is  reflected  as  a  negative  pressure  wave. 

At  a  given  depth  A,  the  time  interval  r(A)  between  the  passage  of  the  front  of  the  shock  wave 
and  the  front  of  the  reflected  wave  is  called  the  “cutoff  time.”  If  for  all  depths  A  the  total 
pressure  of  the  water  immediately  after  cutoff  is  positive,  then  no  cavitation  occurs.  However, 
if  the  total  pressure  is  negative  immediately  after  cutoff,  then  the  water  cavitates.  In  this 
case,  it  is  assumed  that  cavitation  is  present  at  all  depths  A  just  after  cutoff.  When  the  water 
cavitates,  the  water  particles  at  a  depth  A  are  “kicked  off”  with  a  vertical  velocity  11(h) 
which  is  called  the  “kickoff  velocity.”  Using  the  results  presented  in  this  report,  it  can  he 
shown  that  the  cutoff  times  of  interest  are  much  smaller  than  the  times  characterizing  the 
development  of  bulk  cavitation.  Therefore,  it  is  assumed  that  the  time  of  kickoff  is  the  same 
for  all  particles  regardless  of  depth.  A  more  detailed  model  will  be  considered  in  later  work 
on  the  compressible  theory. 

In  the  theory  of  the  present  report,  the  water  is  assumed  to  be  incompressible.  It  is 
also  assumed  that  all  the  water  has  zero  pressure  just  after  cutoff;  this  assumption  is  justi¬ 
fied  because  the  pressure  in  the  cavitated  zone  is  approximately  the  water  vapor  pressure 
which  is  negligible  for  cases  of  practical  interest.  Immediately  after  kickoff,  air  pressure 
decelerates  the  water  particles  at  the  surface.  The  particles  just  below  the  surface  then 
join  up  with  the  surface  particles.  This  is  the  beginning  of  the  formation  of  the  “surface 
slab.”  Air  pressure  continues  to  decelerate  the  surface  slab  so  that  even  more  particles 
below  the  surface  join  the  surface  slab.  In  this  manner,  the  thickness  of  the  surface  slab 
grows  from  zero  initial  thickness.  Figure  1  shows  a  surface  slab  which  has  grown  to  a 
thickness  y  at  the  time  t. 

As  explained  later,  there  is  a  “bottom  zone”  (see  Figure  1)  in  which  the  water 
motionless.  At  time  t,  this  zone  will  have  the  depth  yb.  Since  yb  decreases  as  a  function  of 
time,  this  bottom  zone  moves  upward.  Because  the  thickness  of  the  surface  slab  increases 
as  a  function  of  time,  there  will  be  a  certain  time,  called  the  “closure  time,”  at  which  the 
surface  slab  meets  the  bottom  zone.  After  closure,  the  water  is  assumed  to  have  no  further 
movement. 

KICKOFF  VELOCITY 

The  impulse  imparted  by  a  shock  wave  to  a  column  of  water  of  small  cross-sectional 

area  A  above  a  depth  A  is  /r(A)  AP(T)dT,  where  P(T)  is  the  pressure  of  the  incident  shock 

0 

wave  at  a  time  T  after  arrival.  This  impulse  must  be  equal  to  the  momentum  of  the  column  of 

water  which  is  fh  p All  (h")  dh',  where  p  is  the  density  of  the  water.  Equaling  these  two 
o 

expressions,  differentiating  both  sides  with  respect  to  A,  and  dividing  both  sides  by  pA  gives 
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(1) 


«/(*) 


P(r(h)) 
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drjk) 

dk 


SURFACE  DISPLACEMENT 

Since  the  particle*  in  the  enviuted  tone  are  all  at  taro  pressure,  the  only  deceleration 
they  experience  is  gravity.  Thus,  at  any  time  {'(where  t  *  {'*()>.  the  displacement  of  a 
particle  which  was  kicked  off  at  a  depth  A  is 


(/(A ){'-  j  l'a 

where  g  is  the  deceleration  due  to  gravity.  That  is,  "free  fall"  displacement  will  occur  until 
the  particle  joins  the  surface  slab.  At  the  instant  this  particle  joins  the  surface  slab,  it  will 
have  the  displacement 

t/(A){-  |  {» 

Since  the  water  is  incompressible,  the  thickness  of  the  surface  slab  y  must  be  equal  to  A. 
Thus  the  displacement  of  the  surface  is 


a-t/(A){-  j  {* 
.(/(y){-|  {* 


(8) 


Differentiation  of  Equation  (2]  with  respect  to  (  yields  the  velocity  of  the  surface 


V  - 


dt 


d[V(y)(\ 

dt 


-ft 


■  U(f)  ♦  1 


dV(V) 

df 


(S) 

(4) 
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MOMENTUM  EQUATION 

Now  consider  the  water  in  the  surface  slah  at  time  t  having  a  thickness  y.  At  kick¬ 
off,  all  of  this  water  was  undisplact-d  nnd  had  the  kickoff  velocity  distribution  i  (y*)  for 
y  -  y'-0-  The  momentum  of  n  column  of  small  cross-sectional  area  A  of  this  water  was 

fypAL  (y‘)  dy' 
o 

However,  in  the  time  interval  from  kickoff  to  the  tune  f,  the  momentum  of  the  column  of  water 
of  cross  sectional  area  y|  has  been  reduced  due  to  air  pressure  /*  and  t  gravity  by  the  amount 

raAt  ♦p  Ayyt 

If  there  is  no  momentum  flux  into  (or  out  of)  the  sides  of  this  water  column,  the  momentum  of 
the  column  at  time  t  is 


I y  pAl  (y ')  dy *  -  /»  At  -  pAyyt 
o 

which  must  be  equal  to  pAVy.  That  is, 

pAVy  -  (ypAl'(y')dy'-raAt  -  pAgyt  (51 

o 

which  is  the  momentum  equation.  The  effect  of  momentum  flux  into  (or  out  of)  the  sides  of 
the  water  column  will  be  considered  in  Inter  work  on  the  compressible  theory. 


SURFACE  SLAB  THICKNESS 

Substituting  Equation  (.1]  into  Equation  (5],  dividing  both  sides  by  A,  and  adding 
pyyt  to  both  sides  gives 


i>y 


d\l(y)  ^1 

iit 


pfyU(y')dy’-l>at 

o 


16] 


It  should  be  notea  that  y  is  independent  of  y  in  Equation  |6|.  This  is  expected  since  the  rate 
at  which  the  surface  slab  is  being  decelerated  by  gravity  is  the  same  as  the  rate  at  which  the 
particles  in  the  cavitated  zone  are  being  decelerated  by  gravity. 


pressure  is  negligible  for  eases  of  practical  interest.  In  these  cases,  the  small  error  produced  by  neglecting  the 
water  vapor  pressure  tends  to  be  counteracted  by  the  small  error  caused  by  neglecting  the  increase  in  air  pressure 
due  to  the  sudden  upward  movement  of  the  water  surface. 
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Equation  [oj  can  be  written  as 


or 


dV  (v)  dy 

yt  ^  -  Pi (y')dy'-yl(y)  17) 

ay  dt  o 

fyHy')dy'-yHy)  ~(Pm/p)t 
dy  o 

dt  dV(y)  181 

V  — T- 


It  will  be  shown  later  that  for  all  cases  of  interest  y  is  a  monotonic  increasing  function  of  t. 
Therefore,  we  can  let  t  become  the  dependent  variable  and  y  become  the  independent  variable. 
So,  reciprocation  of  both  sides  of  Equation  (8]  gives 


,  dVM 

dt  "  — 

dy  %V{y')dy'-yHy)-(Pm/p)t 


[9) 


Equation  (0]  relates  the  time  to  the  surface  slab  thickness  and  consequently  governs  the 
motion  of  the  water  surface.  It  can  be  easily  verified  by  direct  substitution  that  the  solution 
of  Equation  [9]  for  which  t  *  0  at  y  -  0  is 


f  *  2  l/rf/(yVy'-yt/(y) )  [10] 

pa  o 

Differentiating  Equation  (10]  with  respect  to  y  and  substituting  the  resulting  expression  into 
Equation  |4]  yields 


t 

V  -  V(y)  -  - - yt 

2  p  y 


[ill 


PHYSICAL  CONSIDERATIONS 

Equation  [  1]  demonstrates  that  for  most  cases  of  interest  U(y)  is  a  monotonic  de- 
creasing  function  of  y.  Thus,  using  Equation  ( 10],  it  can  be  shown  that  dt/dy  ^  0.  That  is, 
y  is  a  monotonic  increasing  function  of  t.  If  (J(y)  were  not  a  monotonic  decreasing  function 
of  y,  several  slabs  would  form  (see  Appendix).  Also,  at  y  ■  0,  dt/dy  ■»  0,  i.e.,  initially  the 
surface  grows  infinitely  fast.  This  can  be  seen  physically  by  remembering  that  initially  the 
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t 


surface  slab  is  of  zero  thickness  and  is  being  infinitely  decelerated  by  the  air  pressure. 

Thus,  the  particles  of  water  just  under  the  surface  will  join  the  surface  slab  immediately. 

Equation  (1]  also  shows  that  as  y  -•  «,  V(y)  +  0  and  fr  U(y*)  dy'  +  constant  for  all 

o 

cases  of  interest.  Therefore,  Equation  ( 10]  implies  that  a s  y  assymptotically, 

i.e.,  as  t  -  tm,  the  slab  grows  infinitely  fast  and  gets  infinitely  thick. 

It  is  now  convenient  to  define  the  surface  displacement  without  the  effect  of  gravity 
as  (see  Equation  (21) 


- V(y)t 


(121 


and  the  surface  velocity  without  the  effect  of  gravity  as  (see  Equation  (11]) 


■  V  ♦  gt 


t 

-  U(V)  -  —  - 
2  P  y 


(IS] 


It  ia  seen  from  the  above  discussion  that  as  y-»«,  $g-»0.  So  it  is  concluded  (Equation  [12]) 
that  initially  rises  from  zero  and  eventually  falls  back  to  zeip  at  t~tm.  Physically,  the 
air  pressure  eventually  reverses  the  upward  motion  of  the  surface  and  forces  it  back  down  to 
its  original  level.  The  fact  that  the  final  surface  level  is  the  same  as  the  initial  surface  level 
is  consistent  with  the  assumption  of  incompressibility.  This  implies  (see  Equation  [13])  that 
Vf  must  decrease  from  its  initial  value  U( 0)  to  zero  and  become  negative.  Finally,  at  t 
Equation  [11]  shows  that  Vg  must  be  zero.  This  means  that  the  acceleration  of  the  surface 
must  be  zero  some  time  after  the  velocity  becomes  negative  but  before 

Because  the  surface  slab  is  assumed  to  be  incompressible  and  of  uniform  velocity,  the 
pressure  in  tho  slab  must  vary  linearly  with  depth  from  P0  at  the  surface  to  PL  just  above  the 
cavitated  zone.  Thus,  the  total  force  per  unit  area  on  the  slab  is  PL  -  Pa  -  pyg  which  must 
be  equal  to  py  dV/dt,  i.e., 


dV 

py  -  pL-pa  -  pyy 


[14] 
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which  becomes 


?L 


PJ_  _/ _ 

4  P  ,  dVM 
*2  — 7— 


(15) 


by  using  Equations  ( 10]  and  [ll].t 

Expanding  V(y)  in  a  Taylor  aeries  about  y« 0  and  aubatituting  Equation  (10)  into 
Equation  (15),  it  can  be  shown  that  as  y-»Q,  P>  +Pa  '4.  It  can  also  be  shown  that  PL  as 
y-*m.  Equation  ( 15]  demonstrates  that  P t  is  always  positive.  If  this  were  not  the  case,  the 
surface  slab  itself  would  be  a  cavitated  cone  which  is  physically  meaningless. 

BOTTOM  ZONE 

It  must  also  be  assumed  that  if  a  particle  which  was  kicked  off  at  a  depth  yb  falls 
back  to  the  same  depth  before  it  joins  the  surface  slab,  then  it  has  no  further  movement. 

That  is,  it  remains  stationary  after  it  returns  to  its  initial  position  rather  than  continuing  to 
fall.  If  this  were  not  the  case,  the  final  water  level  would  be  different  than  the  initial  water 
level.  This  assumption  is  consistent  with  the  assumption  of  incompressibility.  The  time  t 
at  which  this  particle  returns  to  its  initial  depth  is 

2  V(yb) 

/  -  -  lie] 

9 

From  this  equation  and  the  fact  that  V(yb)  is  a  monotonic  decreasing  function  of  yb  and  that 
f/(y6)-«0  as  y6-»»  for  all  cases  of  interest,  it  is  realized  that  at  any  time  t  there  is  a  depth 
below-  which  the  water  is  motionless.  This  deep  cone  of  motionless  water  is  called  "the  bottom 
zone"  (see  Figure  1).  Initially,  at  <- 0,  yb  -  «.  Later,  as  t  increases,  yb  decreases. 

CLOSURE 

Eventually,  there  will  be  a  time  te  at  which  yb  -  y  e  yg.  That  is,  te  is  the  time  at 
which  the  surface  slab  meets  the  bottom  zone.  After  this  time,  all  the  water  is  assumed  to 
be  motionless.  The  time  te  is  called  “the  closure  time,”  and  the  depth  ye  is  called  “the 
closure  depth.”  Substitution  of  Equation  [16]  into  Equation  [2],  when  ybmy*ye%  shows  that 
the  displacement  of  the  surface  at  the  closure  time  is  zero.  This  is  consistent  with  the 


^Alternatively,  can  be  found  by  realizinc  that  it  ia  the  momentum  with  reapect  to  the  aurface  alab  which  ia 
imparted  to  a  unit  area  in  a  unit  time,  i.e.,  =  p  dy/dt  (V (y)  -  V  ),  which  can  be  ahown  to  be  the  a«ne  aa 

Equation  [is]. 
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assumption  of  incompressibility.  Equating  Equations  1 10)  and  1 161,  setting  yfc-y-yr,  and 
multiplying  both  sides  by  y  '*2  gives 


up  \ 

Y -  1  f  e  6(y')dy'-yft(yf)l  -t’(y,)  117) 

*(i  ® 

which  can  be  used  for  determining  ye  and  thus  tf.  Equation  (17)  gives  the  same  values  of 
ye  (and  as  the  model  in  Reference  1. 

Since  L'(yf)  "  yfc'*2,  it  is  seen  by  using  Equation  111)  that  the  velocity  of  the  surface 
slab  at  the  closure  time  is  -(</  4  Pa  \pyr))tr  '2  which  is  also  the  same  as  that  predicl<>d  l»\ 
the  model  in  Reference  |.*  xhis  completes  the  general  exact  solution  of  the  bulk  caviution 
problem  for  this  theory. 

EXPONENTIAL  SHOCK  WAVE 

The  most  important  special  case  of  a  shock  wave  exponentially  decreasing  with  time 
will  now  be  examined.  In  this  case,  the  pressure  history  has  the  form 

P(T)  "  l‘Qr~T  6  lIHl 

where  PQ  is  the  peak  pressure  and  0  is  the  exponential  time  decay  constant.  For  a  plane 
wave  making  an  angle  0  with  the  surface,  the  surface  cutoff  time  is 


r(y)  -  —  y  cos  /3 

c 


19 1 


where  c  is  the  velocity  of  the  shock  wave.  Substituting  Equations  (la)  and  |19|  into  Equation 
(1)  yields 


t  (y)  -  -  cos  fi  r 

,,r 


•  <o  '■"* 


-(7^) 


where 


•_>/» 

'•  77  ™s 


I ‘201 

1211 

1221 


'A  rough  estimate  of  the  closure  pulse  can  be  obtained  by  assuming  that  the  problem  is  strictly  one-dimensional. 
The  closure  pulse  is  then  a  square  wave  which  arrives  al  a  point  of  depth  :  at  the  time  1^4- |z  ~>c|/c,  where  c  is 

the  velocity  of  sound  in  water.  The  duration  of  the  pulse  is  2z/f  for  0  )'e  and  2Xe/c  for  :'>ye.  In  this 

approximation,  z  should  not  be  too  deep.  Perhaps  it  should  be  less  than  about  3 ye.  The  pressure  of  the  pulse  is 
pcf  2  times  the  absolute  value  of  the  velocity  of  the  surface  slab  at  the  closure  time. 
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is  the  kickoff  velocity  of  the  water  surface  and 


k  h 


c$ 


2  cos  fi 


[23] 


is  the  exponential  depth  decay  constant. 

NONDIMENSIONAL  VARIABLES 

The  nondimensional  time  t *  is  defined  as 

and  the  nondimensional  thickness  y*  is  defined  as 

y*  s  y/k 

Using  these  definitions  Equation  [10]  becomes 

f  -  2  [I  -(1  + 

for  the  exponential  wave.  Similiarly,  defining  the  nondimensional  velocity  V*  as 

V *  *  V/UQ 

and  the  nondimensional  gravitational  deceleration  as 

r  ■  9pk/Pa 

Equation  [11]  becomes 


[24] 


[25] 


[26] 


[27] 


[28] 


V *  .  yy*  -  1 _ a*  t* 


2  y* 


[29] 


for  the  exponential  wave.  Finally,  Equation  [2]  determines  the  nondimensional  displacement 


s‘‘s/{k  '«*) 


=  e-Vf  -  —  t* 2 
2 


[30] 


for  the  exponential  wave. 
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CLOSURE  DEPTH  AND  TIME 


Using  Equation  [17],  it  is  found  in  the  case  of  the  exponential  wave 

e>e  =  1  +  y*  +  1  /g*  [31] 

where  y*  s  ye/k  is  the  nondimensional  closure  depth. 

Also,  the  nondimensional  closure  time  is 

C " <<».) '  v) 

■  '(s;)  W2i 


2/y 

l  +  y*  +  l/y* 


Equations  [3 1]  and  [32]  for  the  closure  depth  and  time  are  the  same  as  the  results  of 
of  tho  model  presented  in  Reference  1  as  expected.  This  means  that  the  simple  model  of 
Reference  1  is  just  as  good  as  this  more  detailed  theory  in  finding  the  closure  depth  and 
time.  However,  the  surface  slab  thickness,  velocity,  and  displacement  histones  are  found 
more  accurately  by  using  Equations  [24]  through  130]  of  this  report.  The  experimental  veloci¬ 
ty  and  displacement  histories  in  Reference  l  seem  to  verify  this  more  detailod  theory.  Un¬ 
fortunately  the  accuracy  of  these  histories  is  not  good  enough  to  conclude  that  this  theory  is 
significantly  better  for  predictions  than  the  simplified  model  in  Reference  1. 

CALCULATIONAL  AIDS 

As  an  aid  to  calculations.  Equation  [26]  relating  t*  and  y*  has  been  plotted  in  Figure 
2.  As  explained  before,  y*  grows  infinitely  fast  at  t *  ().  At  t*  0.5285,  the  curve  has  an 

inflection  point.  At  t  =  2,  the  surface  slab  is  again  growing  infinitely  fast  and  is  of  infinite 
thickness. 

Figure  3  shows  the  nondimensional  velocity  and  di -placement  without  the  effect  of 
gravity  (i.e.,  V*  -  F*  -»  y*  t*  and  b*  b*  *(jj*  clotted  using  Equations  [20]  and  [30]) 

as  functions  of  t*.  At  t*  (),  the  deceleration  of  th<  *urfnee  is  infinite.  Later,  at  t*  =  0.715, 
F*  is  zero  and  b*  has  its  maximum  value.  0.201.  Still  later,  at  /*  -  F*  has  its  minimum 

n  A-  p 

value,  -0.218.  Finally,  at  t  2.  and  b ^  ar  •  zero  and  the  acceleration  is  infinite.  Fig¬ 
ure  4  shows  un  enlarged  view  of  the  early  part  of  Figure  3. 
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An  example  of  the  use  of  these  graphs  is  illustrated  in  Figure  5.  In  this  figure,  the 
velocity  history  resulting  from  the  detonation  of  a  10,000-lb  charge  of  HBX-1  at  a  depth  of 
100  ft  and  horizontal  standoff  of  400  ft  is  presented  according  to  the  theory  of  this  report. 

Also  shown  are  the  predictions  in  accordance  with  the  theory  of  Walker  and  Gordon1  and  the 
experimentally  observed  velocity  history.1  The  curves  are  terminated  at  the  time  of  arrival 
of  the  very  complicated  velocity  associated  with  the  expansion  of  the  explosion  bubble  in 
this  particular  experiment.  It  can  be  seen  that  the  theory  of  this  report  gives  a  more  realistic 
characterization  of  the  velocity  history  than  does  the  theory  of  Walker  and  Gordon.1  However, 
there  are  other  experimental  velocity  histories  in  Reference  1  to  which  the  theory  of  Walker 
and  Gordon  appears  to  give  a  better  fit  than  the  theory  of  this  report  at  later  times.  The 
accuracy  of  the  experimental  data  is  not  sufficient  to  determine  that  the  theory  of  this  report 
is  significantly  better. 


SUMMARY  AND  CONCLUSIONS 

A  hulk  cavitation  theory  which  has  a  simple  exact  solution  has  been  presented.  The 
solution  is  physically  plausible  and  is  in  reasonable  agreement  with  the  experimental  data  of 
Walker  and  Gordon.1  Calculations  can  now  be  made  for  the  complete  surface  slab  thickness, 
velocity,  and  displacement  histories  using  Equations  [  10],  [11],  and  [2]  for  any  monotonic 
decreasing  kickoff  velocity  distribution.  The  most  important  special  case  of  the  exponential 
shock  wave  can  be  treated  with  the  aid  of  Figures  2,  3,  and  4.  All  calculations  can  be  per¬ 
formed  without  the  aid  of  a  computer. 
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APPENDIX 

U(y)  FUNCTIONS  WHICH  ARE  NOT  MONOTOWiC  DECREASING 

If  11  (y)  were  not  monotonic  decreasing,  a  slab  other  than  the  surface  slab  would  in* 
slantly  form  starting  at  a  depth  y t  shallower  than  the  point  at  which  6'(y)  begins  to  increase 
and  extending  to  some  depth  y  deeper  than  the  point  at  which  L(y)  is  again  decreasing  (see 
Figure  6).  Since  the  water  is  incompressible,  this  slab  must  be  of  uniform  velocity  and  of 
zero  pressure.  Because  it  must  ha\e  the  same  momentum  as  the  water  in  it  had  initially,  the 

1  r  2 

velocity  of  this  slab  is  —  —  I  i  (y  ')  dy *  which  must  lie  equal  to  l! (y  )  and  to  t’(y,). 

»2-»l  J, 

’I 

The  bulk  cavitation  problem  in  this  case  can  now  be  treated  using  the  procedure  in  the  main 

>  i 

1  r  2 

text  by  replacing  l  (y)  by  -  L'(y')  dy '  -  V  (y.)  -  l  (y,)  for  y  <  y  y,.  When  the 

i  Jv 

•  i 

depth  of  the  surface  slab  is  y it  will  instantly  grow  to  the  depth  y2>  At  this  time,  there  will 
be  a  sudden  increase  in  the  vertical  velocity  of  the  surface.  If  there  is  more  than  one  such 
increase  in  the  U  (y)  function,  these  addi'.onal  increases  can  be  treated  in  a  similar  fashion. 
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SIOMAL  VELOCITY  (WITHOUT  GRAVITY) 


I 


NON  DIMENSIONAL  TIME,  '* 

Figure  4  —  Early  Nondimension  al  Velocity  and  Displacement  without  the  Filled  of 
Gravity  versus  the  Nondimensional  Time 
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NONDIMENSIONAL  DISPLACEMENT  (WITHOUT  GRAVITY) 


Figure  5  -  Comparison  of  Experimental  and  Theoretical  Velocity  Histories 

Charge;  10,000  lb  HBX-1,  burat  depth;  100  ft,  horizontal  atandoff;  400  ft. 


Figure  8  -  A  £/(y)  Function  Which  Is  Not  Monotonic  Decreasing 
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